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Hawking evaporation of Klein-Gordon and Dirac particles in a non-stationary 
Kerr-Newman space-time is investigated by using a method of generalized tor- 
toise coordinate transformation. The location and the temperature of the event 
horizon of a non-stationary Kerr-Newman black hole are derived. It is shown 
that the temperature and the shape of the event horizon depend not only on 
the time but also on the angle. However, the Fermionic spectrum of Dirac par- 
ticles displays a new spin-rotation coupling effect which is absent from that of 
Bosonic distribution of scalar particles. The character of this effect is its obvious 
dependence on different helicity states of particles spin-1/2. 
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I. INTRODUCTION 

The fourth quarter of last century has witnessed various remarkable progress on several 
researches on black hole physics since Hawking's remarkable discovery [|T|. One of these 
aspects is to reveal the thermal properties of many kinds of black holes by miscellaneous 
methods ( [0,^ and references therein). Much efforts have been devoted to the thermal 
radiation of scalar and Dirac fields in some static or stationary black holes In the case 

of a non-stationary axisymmetric black hole, many works have been done on the radiation 
of scalar particles . Though Hawking effect of Dirac particles has been investigated in 
the non-static case [|| , it is difficult to deal with the Hawking evaporation of Dirac particles 
in a non-stationary axisymmetric black hole. The difficulty lies in the non-separability of 
the radial and angular variables for Chandrasekhar-Dirac equations |^ in the non-stationary 
axisymmetry space-time. Recently this dilemma has been attacked by us [0] through con- 
sidering simultaneously the asymptotic behaviors of the first-order and second-order forms 
of Dirac equation near the event horizon. A new term representing the interaction between 
the spin of Dirac particles and the angular momentum of evaporating Kerr black holes was 
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observed in the thermal radiation spectrum of Dirac particles. The character of this spin- 
rotation coupling effect is its obvious dependence on different helicity states of particles 
with spin-1/2. This effect disappears when the space-time degenerates to a spherically 
symmetric black hole of Vaidya-type. It should be noted that this term displayed in the 
Fermi-Dirac spectrum is absent in the Bose-Einstein distribution of Klein-Gordon particles. 

In this paper, we extend the method developed in Ref. |T0[ to tackle with the thermal 
radiation of Klein-Gordon and Dirac particles in a non-stationary Kerr- Newman space-time. 
This analysis is of theoretical interest under current consideration. It is shown that the 
location and the temperature of the event horizon depend on the time and the angle. The 
Fermionic spectrum of Dirac particles displays a spin-rotation coupling effect due to the 
interaction between the particles with spin-1/2 and the black holes with rotation. 

The paper is outlined as follows: In Sec. 2, the location of the event horizon of a non- 
stationary Kerr-Newman black hole is derived by using the method of generalized tortoise 
coordinate transformation. Then Klein-Gordon equation of scalar particles and Dirac equa- 
tion of spinor fields are manipulated in Sec. 3 and Sec. 4, respectively. In Sec. 5, both 
equations for massive particles are recast into a standard wave equation near the event hori- 
zon, and the "surface gravity" of the event horizon is obtained. Sec. 6 is devoted to derive 
the thermal radiation spectra of scalar and spinor particles. In Sec. 7, we present some 
discussions about the spin-rotation coupling effect. 



II. GENERALIZED TORTOISE COORDINATE TRANSFORMATION METHOD 

The metric of a non-stationary Kerr-Newman black hole Il2|,|l3| can be written in the 
advanced Eddington-Finkelstein system as 

, „ A - sin^ ^ + - A . 2 , , ^ , , 

as — — av +2 — asm uavdif — zavar 

Zj Zj 

+2a sin' Odrd^ - J^dO' - + sin^ ^ ^.^2 ^^^2 ^ 

Zj 

where A = r^— 2M(f )r+(5^(t^)+a^, S = r^+a^ cos^ 6 = p*p, p* = r+iacos6, p = r—iacos6, 
and V is the standard advanced time. Both the mass M and the charge Q of the hole depend 
on the time v, but the specific angular momentum a is a constant. The metric (|l]) and the 
electro-magnetic potential 

A=^^(dv-asm'ed^), (2) 

Zj ^ ^ 

is shown by Xu to be an exact solution of the Einstein-Maxwell equations. 



The line element of an evaporating Kerr-Newman black hole is a natural non- 
stationary generalization of the stationary Kerr-Newman solution, but it is of Petrov 
type-II, whereas the latter is of Petrov type-D. The geometry of this space-time is char- 
acterized by three kinds of surfaces of particular interest: the apparent horizons r^j^ = 
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M ± (M^ - Q2 - a^)^/^ the timelike limit surfaces r^^s = M ± (M^ - g2 _ ^2 ^Qg2 ^^,1/2^ 
and the event horizons r^jj = th- The event horizon is necessarily a null-surface r = r(f , 6) 
that satisfies the null-surface conditions g^WiFdjF = and F{v,r,9) = 0. An effective 
method to determine the location and the temperature of the event horizon of a dynamic 
black hole is called the generalized tortoise coordinate transformation (GTCT) which can 
give simultaneously the exact values both of the location and of the temperature of the event 
horizon of a non-stationary black hole. Basically, this method is to reduce Klein-Gordon 
or Dirac equation in a known black hole background to a standard wave equation near the 
event horizon by generalizing the common tortoise-type coordinate r^, = r + {2n)~^ ln(r — r//) 
in a static or stationary space-time |jl5| (where k is the surface gravity of the studied event 
horizon) to a similar form in a non-static or non- stationary space-time and by allowing the 
location of the event horizon rn to be a function of the advanced time v = t + r^, and/or the 
angles 9,{p. 

As the space-time under consideration is symmetric about yj-axis, one can introduce the 
following generalized tortoise coordinate transformation (GTCT) |jlO| 

r, = r + ^ ln[r - rniv, 9)] , 

= V — vq , 9^ = 9 — 9q , (3) 

where = r{v,9) is the location of event horizon, and k, is an adjustable parameter. All 
parameters k, vq and 9q characterize the initial state of the hole and are constant under the 
tortoise transformation. 

Applying the GTCT of Eq. (^ to the null surface equation g^^diFdjF = and then 
taking the r — > rnivo, 9o), v ^ Vq and 9 ^ 9q limits, we arrive at 

- 2{rl + o^Yh + a' sin^ 9,rl + r'^'] (^f)' = , (4) 



in which the vanishing of the coefficient in the square bracket can give the following equation 
to determine the location of the event horizon of an evaporating Kerr-Newman black hole 

/^H - 2{rl + a^)rH + sin^ ^or^ + = , (5) 

where we denote Aj^ = r'jj — 2MrH — Q'^ + a?. The quantities th = dru/dv and r'^ = drH/d9 
depict the change of the event horizon in the advanced time and with the angle, which reflect 
the presence of quantum ergosphere near the event horizon. Eq. (^) means that the location 
of the event horizon is shown as 



M 



M'' + a' sin' 9orjj + r'^^ 2 



2f 



H 



rH = —^± :. ' : ' ^ -a' ' . (6) 



(1 - 2rHy 1 - 2r 



H 



1/2 



The plus (minus) sign corresponds to an outer (inner) event horizon. 
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III. KLEIN-GORDON EQUATION 



In this section, we will consider the asymptotic behavior of minimally electro-magnetic 
coupling Klein-Gordon equation near the event horizon. The explicit form of wave equation 
describing the dynamic behavior of scalar particles with mass /xq and charge q 
1 



'-g 



-.{dk + iqAk 



gg''\d,+iqA,)^ + = ^ 



in the above space-time ([^) is 



92 22^^ 
dr"^ ~^ ^ ~^ ^ dvdr ~^ drdip 

1 92 o . o . 92 



2a-^ — 



sin 9 d(p' 



+ a' sm' 6— + cot 9— + 2r— 

(7f^ OO OV 



d 



+ 2(r - M + iqQr)^ + iqQ - /igS 



$ = 0. 



Under the GTCT 



it can be transformed into 



r^(l - 2rH) - M 



2 Ah - 2rH{rl + 0? 



$ - 2r'rr 

drl dr^dO^ 



+ 2a{l-rH] 



92 



$ + 2(r^ + - f^a^ sin^ 



92 



dr^dv. 
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-2rHri^ - 2iqQrH + r'rr cot 6*0 + + rna sin^ 6'o)— $ = . 



(7) 



(8) 



(9) 



In deriving Eq. (P), we have made use of the event horizon equation to deal with an 
infinite form of 0/0-type to obtain a finite value 

A - 2(r2 + a'^)rH + sin^ 9rj^ + r^^ 



lim 



2{rH - M) - Arur-H 



(10) 



r — Tfj 

By adjusting the parameter k, Eq. (||) can be reduced to a standard wave equation near the 
event horizon. However, we leave it to section 5. 



IV. DIRAC EQUATION 



To write out the explicit form of Dirac equation in the Newman- Penrose (NP) |Tf 
formalism, we establish the following complex null-tetrad system that satisfies the orthogonal 
conditions I n = —m ■ m = 1 

I = dv — a sin^ 6d(p , n = — (^dv — a sin^ Odip'j — dr , 
— I i sin 9[adv — (r^ + a'^)dip] — Srf^j , 

( — i sin 6 adv — (r^ + a'^)d(f —Ild6\. 
V2p L ^ ^ 



m 



m 



(11) 
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and obtain the corresponding directional derivatives 

^ d ^ r'^ + d Ad ad 

1 /. . „ 9 d id 
= — ta sm U— — h T-r H — r 



\^p* ^ dv do sin dip 

S = ^(-iasme— + — -—- — ). (12) 
y^p^ dv do sme dip^ ^ ' 

It is of no difficulty to calculate the non-vanishing NP spin coefficients of the non-stationary 
Kerr-Newman space-time in the above null-tetrad as follows 

1 ia cos 6 rA r — M 

P=—.^ e = 7 



A iasin^ _ {Mr — QQ)iasm.9 

^~2V"' V2p*2 ' V2Sp ' 

_ _ ia sin 6 cot 6 ira sin 9 , . 

a = 71 — p , TT = ^ — , p = — ^ ^ ^ . (13) 

Inserting for the null-tetrad components of electro-magnetic potential 

Qr 

An=—, A ■ I = A ■ m = -A ■ m = , (14) 

Zj 

and the needed spin coefficients into the four coupled Chandrasekhar-Dirac equations in 
the Newman-Penrose formalism 

{D + e-p + tqA-l)Fi + (6 + n-a + iqA- m)F2 = , 

v2 

{A + jj-^f + iqA- n)F2 + {5 + 13-T + iqA- m)Fi = , 

v2 

{D + e* -p* + iqA ■ l)G2 -{6 + r -a* + iqA ■ m)Gi = , 

V2 

(A + /X* - 7* + iqA ■ n)G^ - {6 + (3* - r* + iqA ■ m)G2 = , (15) 

v2 

where /xq, q are the mass and charge of the Dirac particles, respectively, we can get 

/ d r \ 1 / ^ ira sin 9^ 

A / ia cos 9 \ 1 / , a^sin0cos0^ 

)f^.H--^(£t 

d 1 / ^+ zra sin 6'n 



)F2 = 








)Gr = 




)G2 = 
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here we have defined operators 



V=— + A-^\r-M + 2iqQr + 2a— + 2{r^ + a^) — 
or L dip ov 

^ d 1 „ id . . n 9 

^= 7m + o^°^^^^^A lasmO—, 

oO 2 sm 6 dip ov 

1 ^ id ■ . n 9 
0= — + - cote + + lasmO— . 

do 2 sm 6 dip dv 



By substituting Fi = ^Pi, = ^Pa, Gi = ^Qi and G2 = into Eq. (|16D, 

we have 

-^Pi + CP2 = tfioP*Qi , ADP2 + C^Pi = w*Q2 , 

-^Q2-C^Qi = tl^opP2, AVQ,-CQ2 = tiJ.opPi. (17) 

An apparent fact is that the Chandrasekhar-Dirac equations (0) could be satisfied by 
setting 

Qi^P;, Q2^-P*, qQ^-qQ. (18) 



So one may deal with a pair of components Pi, P2 only. Eq. ( \17\ } can not be decoupled 
except in the stationary Kerr-Newman black hole case (M = const) or in the spherical 



symmetry Vaidya-Bonner [|T7l case (a = 0). However, to deal with the problem of Hawking 
radiation, one should be concerned about the asymptotic behavior of Eq. ([T7|) near the 
horizon only. 

First let us consider the limiting form of Eq. (0) near the event horizon. Under the 
transformations (^, Eq. (0) can be reduced to the following forms 

d d 
-{t'h + ia sin ^0^//) jz-Pi + [^h - 2(r^ + a^Yu] TZ-P2 = , 



—Pi + {t'h - la sin e^ru) ^Ps = , (19) 

after being taken limits r — ^ th{vo, 60), v ^ vq and 6 6q. It is interesting to note that a 
similar form holds for Qi,Q2 also. 

If the derivatives ^Pi and ^P2 in Eq. ([ToD are not equal to zero, the existence 
condition of non-trial solutions for Pi and P2 is that the determinant of Eq. (|l^) vanishes, 
which gives exactly the event horizon equation (|^). The relations (|19|) play an important 
role to eliminate the crossing-term of the first-order derivatives in the second-order equation. 
It is consistent to consider the asymptotic behavior of the first-order and second-order Dirac 
equations in the meanwhile because the four-components Dirac spinors should satisfy both 
of them. 

Next we turn to the second-order form of Dirac equations. A direct calculation gives 



6 



d 



(AP— + CC^ - filT^jP^ = /io(asin^Q2 - ^^Ql] 



-ia sin 6 



{2Mr-QQ)^ + M-2iqQr 
or 



P2 



d 



(20) 



Given the GTCT in Eq. @ and after some lengthy calculations, the limiting form of 
Eq. (pOD , when r approaches ruivo, Oq), v goes to vo and 6 goes to 6*0, leads 

+ 2A, - 2Mr, + a)\^ - 2r,^ 



and 



-2a(l - th] 



52 ^2 
+ 2(r|^ + - f-na^ sin^ 6*0) 



{th - 2iqQrH 
d 



— M + ia cos 6'or_H' — ^thTh + cot 6*0 + + r^/a sin Oq) — — >Pi 



-2t{MrH - QQ)a sin ^ot^ A 



-2i(MrjY -gQ)a sin ^0 



+ ia sin (9 
Aj^ - 2{rjj + a2)fj|^ 



Pi 



(21) 



r//(l — 2r 



H) 



M 



+ 2 Ah - 2rH(rjr + a 



^-2r' 



+2a(l - th] 



+ 2{rjj + a- rna^ sin^ 6*, 



- (M - 2iqQrH 



d 



th - ia cos OqTh + r'^ cot + + f/^a sin 6^0) j^-Pj = 



(22) 



where we have replaced the first-order derivative term -^Pi in Eq- ( pl| ) by using the first 
expression of relations (|19|). 



V. HAWKING TEMPERATURE 

In order to reduce Eqs. (|^), (l2l|) and (|^) to a standard form of wave equation near the 
event horizon, we adjust the parameter n such that it satisfies 

rH{l-2rH)-M ^ _ ^.^^^^^ ^ = rj, + a' - r^a' sin^ 9^ , (23) 

which means the "surface gravity" of the horizon is 

rH{l-2fH)-M 



[rjj + a2 - rj^a^ sin^ 6*0) (1 - 2rH) + 2r^ 



2 ' 



(24) 
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where we have used Eq. (^. 

With such a parameter adjustment, these wave equations can be recast into a combined 
form near the event horizon as follows 



92 

h 2 



92 ^ 92 



2C, 



92 



+2(C2 + iCi + iq^u) 



_d_ 



^ = 0, 



where 



ail - th) 



H 



rjj + a'^ — tho? sin^ 6q 

while both Ci and C2 are real, 

-1 



—r 



H 



H 



THd^ sin 6*0 



Co 



2(r|^ + 0.^ — r_f/a^ sin^ ^0 
+rHa'^ sin^ 6'o + cot 6q + 
-1 



r^^(l-4fH)-M + r^ 
2{MrH - QQ)rHa^ sin^ 



2(r?r + — rnO''^ sin^ ^0 



A^,-2rH(r|, + a2) 

2(MrH-QQ)r^a sin^o 



fnCi cos 6*0 — 



(25) 



for = Pi, 

^ _ M -rH + r'lj + rna^ sin^ + cot 
^ 2(r|^ + a2 _ f^a2 gj^2 q^-^ 

^ rua cos 6'p 

' ~ 2{r\ + a2 - f^^a2 sin^ ^p) ' 



for ^ 



P2, and 



C2 



-2rHrH + r'jj cot 6'p + + fj^a^ sin^ 6'p 



2{rl + 0? 



ThO? sin 6^0) 



for \l/ = $. We point out that the above parameter adjustment is a crucial step to achieve a 
standard form of wave equation near the event horizon, which can be viewed as an ordinary 
differential equation because all coefficients in Eq. (pSD are regarded as finite real constants. 
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VI. THERMAL RADIATION SPECTRUM 



Eq. ( p5|) can be treated by separating variables as 

to the radial and angular parts 

B!' + 2(Co + iCi + imVtH + iq^H - ioj)R! = , 9' = AG , (26) 

where A is a real constant introduced in the separation of variables, Cq = AC3 + C2. The 
solutions are 

R = Jl^e2iiu;-,nnH-q^H~C^)r,~2Cor, ^ ^ @ _ ^XB, _ ^27) 

The ingoing wave solution and the outgoing wave solution to Eq. (pSf) , are respectively, 

^ in 5 

Xj/^^^ _ g-«ajD,+jm(/j+A6»*g2i(aj-mrj^-q<I>/f-Ci)r*-2Cor* ^ ^^-^ ^2g^ 

The outgoing wave \E'out is not analytic at the event horizon r = r^, but can be analyti- 
cally continued from the outside of the hole into the inside of the hole by the lower complex 
r-plane 

(r - Th) {th - r)e~''^ 

to 

The relative scattering probability at the event horizon is 



out 



Following the method of Damour-Rufiini-Sannan's |jT5|, the Hawking radiation spectra of 
Klein-Gordon and Dirac particles from the black hole is easily obtained 

1 K 

where m is the azimuthal quantum number, Qh and can be interpreted as the angular 
velocity and electro-magnetic potential of the event horizon of the evaporating Kerr-Newman 
black hole, respectively. In Eq. (0), the upper plus symbol corresponds to the Fermi-Dirac 
distribution, while the lower minus symbol stands for the Bose-Einstein statistics. 
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VII. SPIN-ROTATION COUPLING EFFECT 



The thermal radiation spectra (^) demonstrate that the total interaction energy of 
particles with spin-s in an evaporating Kerr-Newman space-time is 



P ~ ^2 



r'jj + a? — thoP' sin 



mail — th) — pa cos OoTh 



+qQrH + {s + p)MrH-r • 32 

When p = s = 0, it corresponds to the case of scalar fields \E' = in the case of spinor 
fields (s = 1/2), \E' stands for Pi, P2 when p = 1/2, —1/2, respectively. 

The energy spectrum is composed of three parts: Up = uiVLh + q^H + C*!, the first one 
is the rotational energy mVtH arising from the coupling of the orbital angular momentum 
of particles with the rotation of the black hole; the second one is the electro-magnetic 
interaction energy q^H'-, another one is Ci due to the coupling of the intrinsic spin of particles 
and the angular momentum of the black hole, it has no classical correspondence. From the 
explicit expression of the "spin-dependent" term Ci 

Ci = — ^ \-pcos doTH + {s + p)MrH-r "'"^^ff^ J ' (33) 

I-th^ Ah - 2rH{rfj + a^) J 

one can easily find that it vanishes in the case of a stationary Kerr-Newman black hole 
(M = const, th = r'fj = 0) or a Vaidya-type black hole (a = 0, = 0, 7^ 0). 

The term Ci is obviously related to the helicity of particles in different spin states, it 
characterizes a new effect arising from the interaction between the spin of particles and 
the rotation of an evaporating black hole. Because vh and describe the evolution of 
the black hole in the time and the change in the direction, we suggest that the radiative 
mechanism of an evaporating Kerr-Newman black hole can be changed by the quantum 
rotating ergosphere which can be viewed as a mixture of the classical rotating ergosphere 
and quantum ergosphere. 



VIII. CONCLUSIONS 

Equations @ and (^4[) give the location and the temperature of event horizon of a 
non-stationary Kerr-Newman black hole, which depend not only on the advanced time v 



but also on the angle 9. Eq. (|31|) shows the thermal radiant spectra of Klein-Gordon 
and Dirac particles in the non- stationary Kerr-Newman space-time. A difference between 
Bosonic spectrum and Fermionic spectrum appears, that is, a new term Ci in the latter one 
is absent from the former one. The new effect probably arise from the interaction between 
the spin of Dirac particles and the rotation of the evaporating black holes. The feature of 
this spin-rotation coupling effect is its dependence on different helicity states of particles 
with spin-1/2 and its irrelevance to the mass of particles. 



10 



To summarize, we have dealt with Hawking radiation of Klein-Gordon and Dirac particles 
in a non- stationary Kerr-Newman black hole. The spectrum of Dirac particles displays 
another new effect between the spin of the particles and the angular momentum of the hole, 
which is absent from the spectrum of the Klein-Gordon particles. This effect is due to the 
coupling of the intrinsic spin of particles with the rotation of the black holes, it vanishes 
when the space-time becomes a stationary Kerr black hole or a Vaidya-type spherically 
symmetric black hole. This study encompasses previous ones [|T^] (when Q = 0) and |n 
(when a = 0) as special cases. 
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